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of Horizontal Symmetry 

Yasufumi Konishi*) and Ikuo S. SOGAMi**) 
'q^ ■ Physics Department, Kyoto Sangyo University, Kyoto 603-8555, Japan 

^ ; (Received September 28, 2009) 

, We investigate Dirac mass matrices derived in the gauge field theory of a horizontal 

^L(' symmetry generated by a central extension of the Pauli algebra. Through numerical analyses 

D ' of the observed data of the charged fermion masses and the flavor mixing matrix of quarks, 

^0 , values of free parameters in the mass matrices are determined and several empirical relations 

are found among the Yukawa coupling constants. As one specific feature of the theory, we 
find different orderings in squared mass eigenvalues for the up and down quark sectors. 
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§1. Introduction 



The Standard Model (SM) of particle physics possesses no principle to lay any re- 
striction on the pattern of the Yukawa interactions. Nine complex coupling constants 
, are treated as free parameters in every four sectors consisting of three generations 

of quarks and leptons. One way to describe order and variety of the generational 
structure is to postulate a gauge symmetry called the horizontal symmetry In 
^ ! the previous paper, one of the present authors has proposed a gauge field theory 

I of a new horizontal symmetry. The purpose of this paper is to investigate the Dirac 

0^ ■ mass matrices derived in the theory and make numerical analyses of the mass spectra 

i of the charged fermions and the flavor mixing matrix (FMM) of quarks. 

On ' The horizontal (H) symmetry of the theory^) is postulated to be described by 

■ the Lie group generated by a central extension of the Pauli algebra which was found 

i in investigating the FMM of quarks and leptons. ^^'^^ The algebra consisting of four 



generators has the central element identified with the democratic matrix^^ which 
can create hierarchical mass spectra of fundamental fermions. The number of the 
Yukawa coupling constants of the theory is reduced to 4/9 of that of the SM. Through 
' spontaneous breakdown of the symmetry, the theory leads to the Dirac mass matrices 

Ai f for the sector / with definite electric charge. 

The matrix Aif with unique non-Hermitian structure possesses four unknown 
complex parameters. To deduce information on the / sector, we must solve the 
eigenvalue problem of the Hermitian matrix TWjA^j. It should be noticed that 
unspecified quantities involved in these Hermitian matrices are reduced down to ten 
real parameters. Consequently, numerical analysis on the experimental data''^ on the 
six masses and four FMM parameters of the quark sector enables us to determine the 
values of ten real parameters, from which we find several empirical relations among 
the Yukawa coupling constants. 



*' E-mail address: konishi@cc.kyoto-su.ac.jp 
**■* E-mail address; sogami@cc.kyoto-su.ac.jp 
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The Dirac mass matrices deduced from the Lagrangian density of the Yukawa 
interactions are shown expHcitly in §2. We describe formahsms for the eigenvalue 
problem of Ai f^'^f iii §3 and §4. Different mass orderings in the up and down quark 
sectors are explained in §5. In §6, values of the parameters in the mass matrices are 
determined by numerical analysis and specific empirical relations are found among 
the Yukawa coupling constants. Discussion is given in §7 and functional dependence 
among parameters and mass eigenvalues is examined in Appendix. 

§2. Dirac mass matrices 

In the low energy region less than the electroweak scale v = 246 GeV,^) the gauge 
field theory of the H symmetry provides the effective theory for flavor phenomenology. 
Breakdowns of the H and electroweak symmetries lead to the Lagrangian density for 
Dirac masses of the quarks and leptons in the following forms^) 

^M=Y. ^LMf^l+ Yl ^lMf^i + h.c. (2-1) 

f=u,d f=i^,e 

in which ^[ ^ are chiral fermion fields and Mf is the Dirac mass matrix. For the 
up sectors (/ = u, v) of the electroweak symmetry, the mass matrix is given by 

where 



is the democratic element. The four coefficients in the mass matrix are expressed 
in terms of the Yukawa coupling constants Yfi and the vacuum expectation value v 
of the scalar field as a/ = Yfi v, bfi = —1/2 v, 6/2 = v and c/ = 31/4 v. For the 
down sectors (/ = d, e), we find 





Mf = afl + bfi \ +—6/2 1-10 \+CfD (2-4) 



where a/ = Yfi v, 6/1 = Y/2 v, 6/2 = 1/3 ^ ~ 3^/4 

In order to diagonalize the non-Hermitian mass matrix A1/, it is necessary to 
have recourse to the bi-unitary transformation^^ 

^/^-^/^i= -^/diagonal- (2-5) 

To derive the mass eigenvalues and the diagonalizing matrix V^, it is necessary to 
solve the eigenvalue problem for the self-adjoint matrix M. /A1/ as 

MfM\ = m>P'^ (2-6) 
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for each charged fermioii sector (/ = u,d,e). The diagonahzing matrix vf is ob- 
tained in terms of the eigenvectors \v^f^'^) and the FMM of the quark sector is con- 
structed by 



(2-7) 



provided that the eigenvectors are arranged in increasing orders of the masses for up 
and down sectors. 

§3. Eigenvalue problem 1: Quark FMM 

To solve the eigenvalue problem for A4jA4p it turns out convenient to use the 
basis vectors 



1) = :^ 




|2) = 4 




|3) 




(3-1) 



which are eigenvectors of the democratic element D. With these bases, the eigen- 
vector in (2-6) is expanded as 

\v^f^) = xf\l)+yf\2)+zf\3). (3-2) 

For the up quark sector, (2-6) is rewritten for the coefficients of as 
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l^up + 2 


V2C:bui \ 


/ 
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|a„p -\/2a„6*2 




Vu 








Vu 




V 




-V2<6„2 |C„|2 + 2|6„2|' / 


V 




) 
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where Cu = Cu + au- Similarly, for the down quark sector, we obtain 









|adP + 2|6di 



Vd 
\za j 



m 



(d)2 



Vd 



(3-3) 



(3-4) 



for the coefficients oi\v^'^'>), where Cci = Cd + aci + bdi- 

To clarify the counting of independent parameters in these equations and the 
FMM, we define phase factors by 

C*bui = |C„6„i|e^^", a„6;2 = \aubu2K^, aab^^ = |ad6d2|e^^^ Cpdi = \CdbdiK^ 

(3-5) 

and introduce the diagonal phase matrix 



P-^ = diag(e'^/, e^'^f, 1) 



(3-6) 
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for the / sector to adjust phase factors. By separating the diagonal phase matrices, 
the eigenvectors of (3-3) and (3-4) are found, respectively, in the forms 



(u)2 



V2\Cubui\[m^ 
-V^\aubu2\{mf'^'^ - \a. 



P-2KiP) 



V (mf'-Kn(m}"^^-K|^-2|6„in / 

f V2\adba2\imf'-\ad\''-2\bdi\^) \ 

-V2\Cdbdi\{mf' -\ad\^) 



(u)2 



(3-7) 



and 



\ ( 



m 



id)2 



\ad?){m 



{d)2 



(3-8) 



|ad|'-2|6di|2) J 



where m - are eigenvalues of squared masses and ivj are the normalization con- 
stants. Then, with the orthogonal matrices 

0{ = (ul ul ui) (3-9) 

consisting of the vectors Uj, the FMM for the quark sector is calculated to be 



V = ol^poi 

with the diagonal phase matrix 

P = diag( e'^, 1 ) = P"tp<i^ /x = /x^ - 



J^ = l'd-l^u 



(3-10) 



(3-11) 



This FMM includes unknown parameters of eight real numbers and two phases. 
The secular equations for the eigenvalue problems in (3-3) and (3-4) work to fix six 

real parameters in terms of the mass eigenvalues. Consequently, two real numbers 
and two phases remain unspecified in the FMM for the quark sector. 

§4. Eigenvalue problem 2: Mass spectra of charged fermions 

Since the secular equations for both of the eigenvalue problems in (3-3) and (3-4) 
take the same form, the suffix / is omitted for all quantities in this section and the 
Appendix. In terms of the shifted variable s = — |ap, the secular equation is 
obtained as 



- (|C|2 - \af + 2|6|2) - 2{\abf - 2|6i62|^) s + 4|a6i62|^ = (4-1) 



where 



l^r = l^ir + I62I. (4-2) 

Let us solve this equation by the Cardano method. Introducing the dimensionless 
quantities 



^_ l |Cp-|ap + 2|bp ^_ 2 |ap|bp-2|6i62p 
^ |a6i62|3 ^ |a6i62|3 



(4-3) 
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and changing the variable by ,s = ] 0,61 62 (t + P), we obtain the reduced proper 
equation without the second order term as follows: 

^3 _ 3^p2 ^ g^^ _ _ + 4 = 0. (4-4) 

One solution of this equation is determined ast = t++t- by the sum of two quantities 
t+ and t- which are subject to the relations 

4 = 1 (2P^ + 3PQ - 4 ± i^^) (4-5) 

where 

D = -16P2 - 3P2Q2 _ 24pg _ 4q3 + iq (4.5) 

In the analysis below, it is appropriate to use the polar representation = pe*'^' 
and t- = pe-*l^l in which p and 9 are expressed, in terms of P and Q, as 

= ^^^j^^L. ,4.7) 

Then, the three solutions of (4-4) are derived to be 

ti =ojt+ + uH- =2pcos{\e\ + ^) 

t2 = ojH+ + ujt- = 2pcos{\9\ + ^) (4-8) 
= + i- = 2pcos 9 

where to = exp(i27r/3). 

In this way, we have solved the eigenvalue problems in (3-3) and (3-4) obtaining 
the squared masses as follows: 

ml = lap + i - |a|2 + 2|6|2) [l + 2^1 + 5 cos(|e| + f )] , 

ml = lap + i (|(7p - lap + 2|6p) [l + 2vTH^cos(|e| + f )] , (4-9) 

ml = lap + i (|Cp - lap + 2|6p) [l + 2v^Th^cos^] , 

where the parameter 

|aP|6P-2|M.P 
[|CP- |aP + 2|6p]2 ^^^^^ 

is introduced to simplify the expressions. 

As confirmed below, the magnitude of the angle \9\ must be sufficiently small 
for the mass spectra to have hierarchical structure. Here, it is crucially important 
to note that the squared masses have orderings m\<ml< m| and ml < mf < ml, 
respectively, for 6* > and 9 < 0. 

The squared masses in (4-9) depend on the four real quantities jap, |6ip, I62P 
and |Cp. In the present analysis, we interpret inversely that |6ip, |52p and |Cp are 
functions of the masses and the parameter |ap. Then, all quantities for the quark 
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FMM are determined in terms of the observed quark masses and the independent 
parameter I a|. (See the Appendix). 

In the following sections, we make numerical analysis of the qiiark FMM by using 
experimental values of quark masses as inputs and by adjusting the independent 
parameters |a„| and \ad\- 

§5. Mass orderings of the up and down quark sectors 

For numerical analyses below, we use the quark masses at the energy scale of the 
Z boson, i.e., mz = 91.2 GeV. The values calculated by the renormalization group 
equations^°) are given as foUows:^^) 

ruu = l-27taio MeV, rUc = 0.619±0.084 GeV, nit = 171.7±3.0 GeV, 

(5-1) 

ma = 2.90l:lit MeV, = 55^ MeV, nib = 2.89±0.09 GeV. 

The observed FMM of the quark sector has the prominent feature that the ma- 
trix elements decrease rapidly for each step away from the diagonal. To reproduce 
such characteristics, both of the orthogonal matrices and Of^ in (3-9) must ap- 
proximately be close to the unit matrix. Accordingly, as a step for the FMM analysis, 
it is reasonable to examine numerically which of the solutions characterized by > 
and < in (4-9) and the associated eigenvectors in (3-7) and (3-8) can bring the 
orthogonal matrix nearer to the unit matrix. 

Let us apply the positive-^ solution with the experimental mass values in (5-1) 
to examine the orthogonal matrices. For the down quark sector, it is proved that all 
of the diagonal elements of the matrix can approach to the unit for small value 
of \ad\- Contrastingly, for the up quark sector, some of the diagonal elements of the 
matrix O'f^ is shown to be quite smaller than 1 for any value of \au\- The situation 
reverses completely for the negative-^ solution. Numerical calculations with the 
negative-^ solution show that, while approaches to the unit matrix by adjusting 
\au\, of can not be made close to the unit matrix for any value of |ad|. 

Accordingly, it is necessary to choose the positive- and negative-^ solutions, re- 
spectively, for the down and up quark sectors to reproduce the experimental results 
of the quark FMM. The masses and state vectors of the observed down quark mem- 
bers (d, s, b) must be described by the positive-^ solution with normal ordering as 
follows: 

(d) (d) (d) 

nid = ni\ , nis = ni2 , mi, = ; 

1^^) = \v(^n), \vs) = \v('^'^^), \vb) = 

As for the observed up quark members (u, c, t), their masses and state vectors have 
to be identified with the quantities of the negative-6' solutions with partly-reversed 
ordering as follows: 

(u) (u) (u) 

niu = 1712 , nic = ml , mt = , 

1^;^) = |t;W2), \v^) = l-yHl), \vt) = [v^^^^^. 

In the present theory, it is crucial to accept these interpretations of the solutions of 
the eigenvalue problems in (2-6). 
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§6. Hierarchical structure of the Yukawa coupUng constants 

As shown in the Appendix, the quantities |fe/2| and \Cf\ are expressed in 

terms of the masses and the adjustable parameter \af\'^ in the hierarchical approxi- 

mation (771,3 m\ )• Using these results and accepting the interpretations 

in (5-2) and (5-3), we make numerical analysis of the quark FMM. The best fit to 
the FMM data of the Particle Data Group^) is obtained with the following values of 
the four parameters as 

\au\ = 30.4, \ad\ = 13.2, /x = 0.96, u = 2.32 (6-1) 

which lead to the magnitude of the elements of the quark FMM as 

0.974210 0.225705 0.003595 
\V\ = I 0.225473 0.973343 0.041511 | , (6-2) 
0.008723 0.040746 0.999132 

and the Jarlskog invariant measure^^^ 

J = 3.1 X 10-^ (6-3) 

for the CP violation. The values of all parameters for the best fit found above are 
listed in Table I. 

Table I. Values of ten parameters 



up (jiiark (AkA'j down (jiiark (,\lo\') phases 



|au| 


= 30.4 


|Od| 


= 13.2 


H = 0.96 


\bnl\ 


= 831 




= 92.1 


V = 2.32 


\bn2\ 


= 63800 


\bd2\ 


= 818 




\Cu\ 


= 146000 


\C4 


= 2650 





Results in Table I show clearly that the real parameters satisfy the hierarchical 
orderings ^ l^/iP ^ 1^/2 P ^ l^/P each quark sector. By making more 

careful comparison among them, we find approximate relations 

\bui\ , \bd2\ q1 \bu2\ q2 fa A\ 

— - ~ 1, — - ~ 9 , — — ~ 9 (6-4) 

\bd2\ \bdi\ \bui\ 

which result readily in interesting empirical relations 

^^91 ^^92 (6-5) 

among the four Yukawa coupling constants. 

The largest quantity in Table I is \Cf\ for both quark sectors. Note that the 
original parameter Cf in the mass matrix Mf can appear only through the quantity 
Cf, as Cu = Cu + ttu and Cd = Cd + ad + bdi, in the Hermitian matrix MfM^j^. Due to 
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this feature which works to decrease the number of unknown quantities in M.fM.j, 
we can do nothing but determine \cf\ approximately |c„| « \Cu\ and \cd\ ~ \Cd\- 

Using V = 246 GeV and the data in Table I, we are able to fix the magnitudes of 
the Yukawa coupling constants of the quark sectors as in Table II. In the somehow 
crude approximation for |cj|, the constant \Yf^\ have comparatively large uncertainty. 



Table \1. Yukawa coupliiig const ant s 



/ 


u 




d 


e 




\Yfi\ 


1.2 X 10" 


-4 


5.4 X IQ-^ 


3.5 X 10" 


-5 


\Yf2\ 


3.4 X 10" 


-3 


3.7 X 10""* 


4.1 X 10" 


-4 




2.6 X 10" 


-1 


3.3 X lO^'^ 


3.6 X 10" 


-3 


\Yf4\ 


2.0 X 10" 


1 


3.6 X 10^3 


1.6 X 10" 


-3 



Thus far, the Yukawa coupling constants of quark sectors are calculated so as 
to recreate the observed data of the quark FMM. This approach is not applicable, 
as it stands, to the lepton sector. Here, let us find an approximate scheme with less 

numbers of adjustable parameters based on the benefit of hindsight on the empirical 
relations in (6-5) and apply it to analyze the masses of charged leptons. 
For its purpose, we introduce a new variable /3 by the equations 

|^di| = 9/3, |&d2| = |6«i| = 9'/3, |6«2| = 9^/3, (6-6) 

which verify all of the relations in (6-4), and reinvestigate the quark masses and FMM 
in terms of seven quantities |a„|, |arf|, |C„|, |Cd|, /i, and 13. It turns out possible to 



Table 111. X'alues of parameters for approximate estimation 



up quark (MeV) 


down quark (MeV) 


common parameters 


\au\ = 30 


fid = 13 


^J. = 0.95, z/ = 2.3 


\Cu\ = 142000 


|Cd| = 2620 


/? = 10 



explain all observed data within the range of experimental errors. In fact, by using 
the values of seven quantities in Table III, we obtain 



ruu = 1.20 MeV, 
nid = 2.87 MeV, 



mt = 169.6 GeV, 
rrib = 2.86 GeV, 



0.628 GeV, 
- 53.8 MeV, 

for the six quark masses and 

/ 0.973984 0.226589 0.003509 \ 
\V\ = 0.226452 0.973161 0.040962 

y 0.008616 0.040199 0.999155 j 

for the magnitudes of the quark FMM elements and the Jarlskog parameter. 



J = 2.9 X 10-^ 



(6-7) 



(6-8) 



Dirac Mass Matrices in Gauge Field Theory of Horizontal Symmetry 



9 



For both sectors of the down quark and the charged lepton, the Dirac mass 
matrices take the common forms in (2-4). To analyze the masses of the charged 
leptons in this approximate scheme, we assume that the empirical equations related 
to the down quarks in (6-4) hold also in the charged lepton sector. Then, following 
the equations in (6-6), we introduce a new variable (3e by the relations 

\bel\=9Pe, \be2\=9^Pe (6-9) 

and analyze the charged lepton masses in terms of the three parameters |ae|, |Ce| 
and Pe- Numerical estimation shows readily that the charged lepton masses^'')' 

me = 0.4866MeV, = 102.7 MeV, = 1746MeV (6-10) 

can be reproduced by adjusting the parameters as follows: 

\ae\ = 8.537, \Ce\ = 1197, /?e = 11.06 (6-11) 

from which the Yukawa coupling constants for the charged lepton sector can be fixed 
provided that |ce| ~ \C(,\. The results are included in Table II. Fig.l shows behavior 
of the Yukawa coupling constants of charged fermion sectors. 




Fig. 1. The Yukawa coupling constants \Yfi\, \Yf2\ and |i/3|. 



§7. Discussion 

We have investigated the non-Hermitian mass matrices Mf in (2-2) and (2-4) 

deduced in the gauge field theory of the horizontal symmetry generated by the cen- 
tral extension of the Pauli algebra. While the matrix Ai / possesses four complex 
unknown parameters, the Hermitian combination TWyA^J- depends on four real num- 
bers and two phases. Owing to this fact, we are able to make numerical analyses on 
the problem of the six quark masses and the quark FMM with four parameters. 

To solve the problem effectively, six real parameters are interpreted as functions 
of quark masses and remaining two parameters. Using the observed masses as in- 
puts and adjusting values of two independent parameters and two phases, we have 
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reproduced the observed FMM and succeeded to calculate the Yukawa coupling con- 
stants. By careful examination of the numerical results, empirical relations in (6-5) 
are found among the coupling constants for the quark sector. It is beyond the scope 
of the present theory to elucidate physical implications of those relations. 

The estimated values of the Yukawa coupling constants in Table II show novel 
orderings with sizable amount of variation. It should be noticed that the largest 
ratio between the coupling constants, |yu3|/|yui| ~ |^^u4|/|i^ui| ~ 10^, is considerably 
smaller than that between the observed quark masses mt/m 10^ This means 
that varieties observed directly in the low energy flavor physics are much reduced at 
the level of the coupling constants of the Yukawa interactions. 

All results of numerical analyses so far were obtained in the hierarchical approxi- 
mation (mg ^ m\ , m2 ) by using the formulae in the Appendix which express 
1^/2! and \Cf \ as functions of |a/p. It should be noticed, however, that we are 
able to find almost the same results by adjusting the values of real quantities 
\bf2\, \Cf \ and |a/p and two phases fi and u. 

The squared masses of the charged fermions are derived in the exact formula 
in (4-9). Depending on the sign of the angle 9, the formula has the solutions with 
normal and partly-reversed orderings in (5-2) and (5-3). As confirmed in §5, we 
have to describe the down and up quark states by the normal and partly-reversed 
solutions, respectively, in order to reproduce the observed FMM. This is one of 
unique features of the present formalism for the horizontal gauge symmetry. 

We made a speculative analysis of the charged lepton masses by postulating 
hypothetical relations (6-9) as analogues of the equations for the down quark sector 
in (6-6). To proceed a full investigation on the charged and neutral lepton sectors, 
however, it is required to solve the combined eigenvalue problem which involves not 
only the Dirac mass matrices but also the Major ana mass matrix deduced for the 
neutrino sector.^) We will study this problem in near future. 

Appendix A 

Dependence of\bi\, I62I and \C\ on \a\ 



In this scheme, we interpret the real quantities |6i|, I62I and |C| as functions of 
the quark masses and the parameter \a\. For its purpose, it is appropriate to define 
a mass M for reference by the relation 

mi + mi + m? = \Cf + 2\bf + 2\af = 3M^ (A-1) 

Evidently observed data in (5-1) show hierarchical orderings of squared masses, 
M^, m| S> m^, mf, for each sector. To deduce convenient relations for numerical 
analyses, we express all quantities as power series of M^. 

In the hierarchical limit, P and |Cp take very large values. Equation (4-7) 
shows that (tan3|^|)^ can be decomposed in the power series of M^^, leading to the 
approximate relation 

(M)^.4^ + lWWz*feQ!. (A.2) 

^ ^ M6 3 ^ ' 
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Similarly, the quantity 6 defined by (4-10) is approximated by 

2 |aP|bP-2|M.P 
•^-3 M4 • ^^^^ 

Decomposition of the eigenvalues and m| in (4-9) with respect to results in 
the following expressions 

ml + mlc^ 2|ap - - ml 2V3M^e. (A-4) 



Eliminating 9 and S from these equations, we obtain the relations 



Then, substitution of the last expression for |6p into the defining equation (A-1) 
allows to express |Cp in the form 



2 2 



|C|2 = 3^^m2 - 2|g|2. (A.6) 



Finally, to distinguish between \hi\ and |52|, we introduce an extra quantity k by 

= - Mk, = \\h\^ + Mk. (A-7) 

From (A- 5), k is calculated to be 



'3 (|a|^ — mlmf)^ 
4|a|2 V4 K 

Consequently, the quantities |C|, |6i| and I62I are determined as the functions 

of quark masses and the parameter \a\ for each sector. The quark FMM is now 
expressible in terms of the observed mass values and the four adjustable parameters 
|au|, |arf|, 11 and v. 
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